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Abstract. In this work, we consider an inverse problem of determining a source 
term for a structural acoustic partial differentia equation (PDE) model, comprised 
of a two or three-dimensional interior acoustic wave equation coupled to a Kirchoff 
plate equation, with the coupling being accomplished across a boundary interface. 
For this PDE system, we obtain the uniqueness and stability estimate for the source 
term from a single measurement of boundary values of the "structure" . The proof 
of uniqueness is based on Carleman estimate. Then, by means of an observability 
inequality and a compactness/uniqueness argument, we can get the stability result. 
Finally, an operator theoretic approach gives us the regularity needed for the initial 
conditions in order to get the desired stability estimate. 
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1. Introduction and Main Results 



1.1. Statement of the Problem. Let Q be an open bounded subset of M? or M 3 
with smooth boundary T of class C 2 , and we designate a nonempty simply connected 
segment of T as T with then r = T U L\ and r fl L^ = 0. We consider here the 
following system comprised of a "coupling" between a wave equation and an elastic 
plate-like equation: 



ztt(x,t) 
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Az(x, t) + q(x)z(x, t) 


-v tt {x, t) - A 2 v(x, t) - A 2 v t {x, t) 

x,t) = 
: v t (x,t) 
z (x) 
zi(x) 
v (x) 
v x {x) 



in VI x [0, T] 
on Ti x [0, T] 
on T x [0, T] 
on dT x [0,T] 
on T x [0,T] 
in Q 
in Q 
on r 
on r 



where the coupling occurs across the boundary interface r . [zq, z±, Vo, Vi] are the 
given initial conditions and q(x) is a time-independent unknown coefficient. For this 
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system, notice that the map {q} — > {z(q),v(q)} is nonlinear, therefore we consider 
the following nonlinear inverse problem: Let {z = z(q),v = v(q)} be the weak 
solution to system (11.11) . Under suitable geometrical conditions on Ti = T \ T , is it 
possible to retrieve q(x), i6 0, from measurement of Vu{q) on To x [0, T]? In other 
words, is it possible to recover the internal wave potential from the observation of 
the acceleration of the elastic plate. 

Our emphasis here that we determine the interior acoustic property from observing 
the acceleration of the elastic wall (portion of the boundary), is not only due to phys- 
ical consideration, but also to the implications of such inverse type analysis related 
to the coupling nature of the structural acoustic flow. In many structural acoustics 
applications, the problem of controlling interior acoustic properties is directly cor- 
related with the problem of controlling structural vibrations since the interior noise 
fields are often generated by the vibrations of an enclosing structure. An important 
example of this is the problem of controlling interior aircraft cabin noise which is 
caused by fuselage vibrations that are induced by the low frequency high magnitude 
exterior noise fields generated by the engines. 

The primary goal in this paper is to study the uniqueness and stability of the inte- 
rior time-independent unknown coefficient q(x) in some appropriate function space. 
More precisely, we consider the follow uniqueness and stability problems: 

Uniqueness in the nonlinear inverse problem 

Let {z = z(q),v = v(q)} be the weak solution to system (II. ip . Under geometrical 
conditions on Ti, does the acceleration of the wall Vtt\r x[a,T\ determine q(x) uniquely? 
In other words, does 



imply q(x) = p(x) in Q7 

Stability in the nonlinear inverse problem 

Let {z(q),v(q)}, {z(p),v(p)} be weak solutions to system (II. ip with corresponding 
coefficients q(x) and p(x). Under geometric conditions on r 1; is it possible to estimate 
Ik - Ph 2 (n) by some suitable norms of {v tt (q) - %(p))|rox[o,T]? 

In order to study the nonlinear inverse problem, we first linearize (ll.ip and hence 
we consider the following system: 



v t t(q)\r x{o,T] = v tt (p)\r x[o,T} 




Aw(x, t) — q(x)w(x, t) = f(x)R(x, t) 


—u u {x, t) — A 2 u(x, t) — A 2 u t {x, t) 


u t {x,t) 



infix [0,T] 
on Ti x [0,T] 
on T x [0,T] 
on dT x [0, T) 
on T x [0,T] 
in fl 



w t (x,t) 
u(x, t) ■ 



du 



(1.2) 
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where q G L°°(0) is given, R(x,t) is fixed suitably while f(x) is an unknown time- 
independent coefficient. For this linearized system, we have the advantage that the 
ma P {/} — >* { w if)i u if)} is linear, hence we consider the corresponding linear inverse 
problem: 

Uniqueness in the linear inverse problem 

Let {w = w(f), u = «(/)} be the weak solution to system (11. 2p . Under geometrical 
conditions on Fx, does Uu\r x[o,T] determine f(x) uniquely? In other words, does 

Utt(f)\r x[o,T] = 

imply f(x) = in f2? 

Stability in the linear inverse problem 

Let {w = w(f),u = u(f)} be the weak solution to system (11. 2p . Under geomet- 
rical conditions on Fx, is it possible to estimate ||/||L 2 (n) by some suitable norms of 

u tt |r x[o,T]? 

Remark 1.1. In our models (11.11) and (11.2j) we regard t = ? as the initial time. This 
is not essential, because the change of independent variables t — > t — — transforms 
t — ~ to t = 0. However, this is convenient for us to apply the Carleman estimate 
established in [29]. In fact, one can keep t = as initial moment by doing an even 
extension of w and u to Q x [— T, T], but then the Carleman estimate in [29] needs 
to be modified accordingly. 

1.2. Literature and Motivation. The PDE system (ll.ip is an example of a struc- 
tural acoustic interaction. It mathematically describes the interaction of a vibrating 
beam/plate in an enclosed acoustic field or chamber. In this situation, the boundary 
segment F\ represents the "hard" walls of the chamber Q, with To being the flexible 
portion of the chamber wall. The flow with in the chamber is assumed to be of acous- 
tic wave type, and hence the presence of the wave equation in Q, satisfied by z in (ll.ip . 
coupled to a structural plate equation (in variable v) on the flexible boundary por- 
tion r . This type of PDE models has long existed in the literature and has been an 
object of intensive experimental and numerical studies at the Nasa Langley Research 
Center pT|[9^[Tu]. Moreover, recent innovations in smart material technology and the 
potential applications of these innovations in control engineering design have greatly 
increased the interest in studying these structural acoustic models. As a result, there 
has been a lot of recent contributions to the literature deal with various topis; e.g., 

optimal control, stability, controllability, regularity [D0|3lll|5j|6l[71|8l[Tl|23]. 
However, to the best of our knowledge, there are no results available in the literature 
for our inverse type analysis on the model. 

On the other hand, the interest to the inverse problem has been stimulated by the 
studies of applied problems such as geophysics, medical imaging, scattering, nonde- 
structive testing and so on. These problems are of the determination of unknown 
coefficients of differential equations which are the functions depending on the point 
of the space [HI [151 US] • F° r the uniqueness in multidimensional inverse problem with 
a single boundary observation, the pioneering paper by Bukhgeim and Klibanov [12] 
provides a methodology based on a type of exponential weighted energy estimate, 
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which is usually referred as the Carleman estimate since the original work [13] by 
Carleman. After [12] . several papers concerning inverse problems by using Carleman 
estimate have been published (e.g. p2J[2l]). In particular, for the inverse hyperbolic 
type problems that is related to our concern in this paper, there has been intensively 
studies [HI [19j |20j [32j [37]. However, we mentioned again that there is not any such 
uniqueness and stability analysis for the structural acoustic models or even in general 
coupled PDE systems. This motivates the work of the present paper. 

The usual problem setting for inverse hyperbolic problem includes determining a 
coefficient from measurements on the whole boundary or part of the boundary, either 
Dirichlet type [H HU El EZ] or Neumann type [TgJ □!]. Usually the coefficient 
describes a physical property of the medium (e.g. the elastic modulus in Hooke's 
law), and the inverse problem is to determine such a property. In our formulation 
of the inverse problem, we need to determine the time-independent wave potential 
q(x) by observing the acceleration from the flexible portion of the boundary Tq. 
The mathematical challenge in this problem stems from the fact that we are dealing 
with the "coupling" on the part of the boundary and the main technical difficulty 
associated with this structure is the lack of the compactness of the resolvent. As a 
result, the space regularity for the solution of the wave equation component is limited 
by the structure on the plate and hence this will prevent us going to higher dimension 
(n > 7) no matter how smooth the initial data is. This is a distinguished feature 
of this structural acoustic model comparing to the purely wave equation model as 
in that case the solution can be as smooth as we want as long as the initial data 
is smooth enough. In this present paper, we prove the cases where the dimension 
n = 2 and 3 (physical meaningful cases) by using the Carleman estimate for the 
Neumann problem in [29J and an operator theoretic formulation. We show that 
indeed the observation of the acceleration on the plate can determine the potential 
q under some restrictions on the initial data and some geometrical conditions on the 
boundary. As we mentioned, the argument will also work for dimension up to n = 7. 

1.3. Main Assumptions and Preliminaries. In this section we state the main 
geometrical assumptions throughout this paper. These assumptions are essential in 
order to establish the Carleman estimate stated in section 2. 

Let v = [fx, ■ ■ • ,v n ] be the unit outward normal vector on T, and let = V ■ 
v denote the corresponding normal derivative Moreover, we assume the following 
geometric conditions on Y x = T \ T : 

(A.l) There exists a strictly convex (real-valued) non-negative function d : O — > 
R + , of class C 3 (Q), such that, if we introduce the (conservative) vector field h(x) = 
[hi(x), ■ ■ ■ , h n (x)} = Vd(x),x G Q, then the following two properties hold true: 



(i) 




dd 
dv 



Vd- v = h ■ v = 0; h = Vd 
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(ii) the (symmetric) Hessian matrix Hd of d(x) [i.e., the Jacobian matrix of h(x)} 
is strictly positive definite on Q: there exists a constant p > such that for all x E Q: 



;i.4) 



H d {x) = Jh(x) 



dxixi 



l d XnXl 



(A. 2) d(x) has no critical point on Q: 
(1.5) inf|%)l = 



dxiX n 

dx^ x™ 



inf \Vd(x) 

xeQ 





- dhi 


dh x - 




dx\ 


8x n 




dh n 


dh n 






dx n _ 


= s 


> 





>pl 



Remark 1.2. One canonical example is that is flat (not the case in our problem 
setting here), where then we can take d(x) = \x — xq\ 2 , with xq on the hyperplane 
containing Ti and outside Q, then h(x) = Vd(x) = 2(x — x Q ) is radial. However, in 
general h(x) is not necessary radial. In particularly in our case where Ti is convex, the 
corresponding required d(x) can also be explicitly constructed. For more examples of 
such function d(x) with different geometries of Ti, we refer to the appendix of 



Next we introduce an abstract operator theoretic formulation associated to (II. ip 
for which we will need the following facts and definitions: Let the operator A be 



'1.6) Az=-Az-q(x)z, D(A) = {z : Az + q(x)z G L 2 (Q), — 

ov 



0} 



Notice the lower-order part is a perturbation which preserves generation of the self- 
adjoint principle part A N (e.g. [2TJ ) , where A N : L 2 (Q) D D(A N ) — > L 2 (Q) is defined 
by: 



(1.7) A N z = -Az, D{A N 
Then A^ is positive self-adjoint and 

(1.8) D{a\) 
Then we define the Neumann map N by: 



H7 

{z : Az G L 2 (£l), — 



0} 



Hy\ x {n) = {z:ze H\n),— = on 



(1.9) 

By elliptic theory 
(1.10) 

Now we define 
(1.11) 



Az = in Q 
Vy { % on Y x 
M = 9 on T 

NeC(L 2 (T ),H 3 r /2 (Q)) 

B = A N N : £ 2 (r ) — )■ d(a\ [ )' 
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via the conjugation B* = N*A N . Then with v G L 2 (T) and for any y G D(A^) we 
have 

(1.12) - (B*y, v) r = ~(N*A N y, v) T = ~(A N y, Nv) n = (Ay, Nv) n 

= (y, A(Nv)) n + (^, Nv) r - (y, ^^) r = -(y, v) Fo 

dy i 
by Green's theorem, the definition of N and the fact — = on r\ when y G D(A 2 N ). 

In other words, we have 

(1.13) N*A N y 



y, 



on r 



o 



o, on r\ 



for y G D(Al) 



i.e. B* = N*An is the restriction of the trace map from if 1 (ft) to (To). 
Last we set A : L 2 (T ) D D(A) ->■ L 2 (r ) to be 

(1.14) A = A 2 , D(A) = {ve H 2 (T ) : A 2 v G L 2 (r )} 

where H$(r ) = {v e H 2 (Q) : d = |^ = on <9r }. A is self-adjoint, positive 
definite, and we have the characterization 



(1.15) 
Now set 

(1.16) 



D(A*) = H 2 (T ) 



A 





-A N + q 





I 



-B* 







-A 




B 
I 

-A 



on the energy space 

17) H = D(4) x ^ 2 (^) x £>(A*) x L 2 (r ) 

= //^(ft) x L 2 (ft) x # 2 (r ) x L 2 (r ) 

Then we have the domain of the operator A 
(1.18) 

£>0A) = {[^o,^i,^o,^i] r e [^(^l)] 2 x [D(A^)] 2 such that 
- z + Nvx G D(A N ) and t> + t>i G D(A)} 
= {[z ,z u v Q , Vl ] T : 2: e ^(ft),^ G H^(Q),v G # 2 (r ),t;i G # 2 (r ), 

(A + g)z G L 2 (ft), ^ = ^011 T and v + vt e D(A)} 

= {[z ,z 1 ,v ,v 1 ] T : z G H 2 i (n),z 1 G H^(tt),v G # 2 (r ),t;i G # 2 (r ), 

= vi on T and w + v± G £>(A)j 

Of 
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where in the last step we get z G H 2 (Q) from q G L°°(fl) and (A + q)z G L 2 (fi) due 
to elliptic theory. Therefore with these notations, the original system (11 .ip becomes 
to the first order abstract differential equation 

(1.19) 1 = A V 

where y = [z, z t , v, Vt] T . From semigroup theory, when the initial conditions [zq, Z\, Vq, V\] 
are in D(A) we have that the solution y satisfies 

(1.20) yeD(A), y t eH 

Remark 1.3. The structure of A reflects the coupled nature of this structural acoustic 
system (II. ip . One distinguished feature of the system is that the resolvent of A is 
not compact. However, it can still be shown that A generates a Co-semigroup of 
contractions {e At } t >o which establishes the well-posedness of the system [3]. 

1.4. Main results. For the inverse problems stated in section 1.1, we have the fol- 
lowing results: 

Theorem 1.4. (Uniqueness for the linear inverse problem) Under the main assump- 
tions (A.l), (A. 2) and let 



(1.21) 

Moreover, let 

(1.22) 

and 

(1.23) 



T > 2^ /maxd(x) 



R G W 6 >°°(Q) 



T 

R ( x, — 
2 



> r > 0, 



T 

Rt ( x, — 



>n>o 



for some positive constants r , r\ and iGO. In addition, let 

(1.24) q G L°°(Q) 

If the weak solution {w = w(f),u = u(f)} to system f ll.2p satisfies 

(1.25) w, w t , w tt G H 2 (Q) = H 2 (0, T'L 2 (n)) n L 2 (0, T; H 2 (n)) 
and 

(1.26) u u (f)(x,t) =0, xeT ,te[0,T] 
then f(x) = 0, x G f2. 

Theorem 1.5. (Uniqueness for the nonlinear inverse problem) Under the main as- 
sumptions (A.l), (A. 2), assume (jl.2ip and 

(1.27) q, P eL°°(n) 
Let either of z(q) and z(p) satisfy 

(1.28) z G W 3 '°°(Q) 
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Moreover, let 

(1.29) \zq{x)\ > s > 0, \ Zl (x)\ > s x >0 

for some positive constants sq, s± and x G Q. If the weak solutions {z(q), v(q)} and 
{z(p),v(p)} to system fll.il) satisfy 

(1.30) z(q) - z(p),zt(q) - z t (p),z tt (q) - z tt (p) G H 2 (Q) 
and 

(1.31) v tt (q)(x,t)=v tt (p)(x,t), xeT ,te[0,T] 
then q{x) = p(x), x G Q. 

Theorem 1.6. (Stability for the linear inverse problem) Under the main assumptions 
(A.l), (A. 2), assume fOTj) . (02]) . ffL23|) and flTMj) . Moreover, let 

(1.32) R t e H^ +e (0,T;L°°(n)) 

for some < e < |. Then there exists a constant C = C(Q,T,T ,(p,q, R) > such 
that 

(1.33) ||/||za(n) < C (\\u tt \\ L 2 {Tox[0iT]) + ||wttt||L2 (rox[0|T]) + ||A 2 Mit|U 2 (r x[o,T])) 
for all f G L?(VL). 

Theorem 1.7. (Stability for the nonlinear inverse problem) Under the main assump- 
tions (A.l), (A. 2), assume ffPTD . (Q7|) . (OS]) and Moreover, Zet t/ie initial 
data satisfy the compatibility condition 

(1) When n = 2, [z , z±, v , vi] G D(A 2 ) where 

D(A 2 ) = {[z ,z 1 ,v ,v 1 f : z G H* 1 (ty,z 1 G ^(O),^ e #0^0)^1 e #0^0), 

i(^ + Wl ) + g # 2 (r ), ^ - i4(v + Vl ) - B* Zl G D(i4), 
i dzi , 

a^' r ° = Vu a^' r ° = ~ ^° + Vl > B Zl * 

(2) VF/ien n = 3, [z , z 1 ,v ,Vi] G _D(*4. 3 ) where 

D(A 3 ) = {[z ,z 1 ,v ,v 1 ] T : z G hI(Q), Zi G H^(Q),v G # 2 (r ),*>i G # 2 (r ), 
A(v + Vl ) + B*z x G # 2 (r ), -rr^lro = vt, — -|r = -;i(vo + «i) - 
A(w + uj) + B*z x + - i(v + vx) - + B*[(-A N + g)^ + Bvi] G 
4(u a - 4(^o + *>i) - 5*^) + S*[(-A,v + g)* + e # 2 (r ), 
g[(z^^+^ll | ro = _ A (v + - 5*,,] - + g), + Bv l]} 

Then there exists a constant C = C(Q, T, T , ip, q,p, z , z\, v , vi) > such 
that 

(1.34) ||g-p||La ( n) < C (\\v tt (q) - v tt (p)\\mr x[o,T}) + \\v ttt (q) - v m (p) \\&(r x[o,T\) 

+ ||A 2 (^(g) - %(p))||L2 (rox[0>T]) ) 
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for all q,p G W 1,0 °(Q) when n = 2 and all q,p G VF 2 '°°(fi) u>/ien n = 3. 

The rest of this paper is organized as follows: In section 2 we give the key Carleman 
estimate that is used in the proof of uniqueness result. Based on the same Carleman 
estimate, we also prove an observability inequality that is needed in section 5. Section 
3 to 6 are devoted to the proofs of our main results Theorems 11.41 to 11.71 Some 
concluding remarks will be given in section 7. 

2. Carleman estimate and observability inequality 

2.1. Carleman Estimate. In this section, we state a Carleman estimate result that 
plays a key role in the proof of our uniqueness theorem. The result is due to [29] . 
We first introduce the pseudo-convex function <p(x, t) defined by 

(2.1) (p(x,t) = d(x)-cU-^\ ; iGO,te[0,T] 

where T is as in (jl.2ip and < c < 1 is selected as follows: By (jl.2ip . there exists 
5 > such that 

(2.2) T 2 > 4 max d{x) + 45 

For this 5 > 0, there exists a constant c, < c < 1, such that 

(2.3) cT 2 > 4 max d(x) +45 

Henceforth, with T and c chosen as described above, this function tp(x, t) has the 
following properties: 

(a) For the constant 5 > fixed in ( 12.21) and for any t > 

(2.4) <p(x,t)<(p(x,^), <p{x,0) = (p{x,T) < d{x) - < -5 

uniformly in x G O. 

(b) There are t and t 1 , with < to < \ < ti < T, such that we have 

(2.5) _min <p(x, t) > a 

xen,te[to,*x] 

where < a < mm x< --^ d(x) . 

Moreover, if we introduce the space Q(o~) that is defined by the following 

(2.6) Q{p) = {{x,t)\x eQ,0<t< T,(p(x,t) > a > 0} 
Then an important property of Q(o~) is that (see [29|): 

(2.7) [t , tjxflc Q(a) C [0, T] x Q. 
Then for the wave equation of the form 

(2.8) w a (x,t) - Aw(x,t) -q(x)w(x,t) = F(x,t), xGfi,tG[0,T] 
we have the following Carleman-type estimate: 
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Theorem 2.1. Under the main assumptions (A.l) and (A. 2), with (f(x,t) defined 
in (12. ip . Let w G H 2 (Q) be a solution of the equation (12. 8p where q G L°°(Q) and 
F G L,2(Q). Then the following one parameter family of estimates hold true, with 
p > 0, (3 > 0, for all r > sufficiently large and e > small: 

(2.9) 

BT\ w + 2 [ e 2Ttfi \F\ 2 dQ + C ltT e 2ra [ w 2 dQ > (rep-2C T ) [ e 2r * (w 2 + \Vw\ 2 ) dQ 
Jq Jq Jq 

+ (2r 3 /3 + 0{t 2 ) - 2C T ) [ e 2rip w 2 dxdt - c T T 3 e - 2TS [E w (0) + E W (T)} 

JQ(ct) 

Here 5 > 0, a > are the constants in (12. 2]) , (12.51) . while Ct, ct and C\ t T are positive 
constants depending on T and d. In addition, the boundary terms BT\ W are given 
explicitly by 



BT\ 



2r 



e 2rLp [w 2 t - \Vw\ 2 )h ■ udTdt 



-To 
T 



+ 8cr 
+ 



T dw . 
— )w t -—di dt 
2 J dv 



(2.10) 



Ur (t - 

4T / J e 2rv (h-Vw)^dTdt 



+ 4r 2 



+ 2r 



\h\ 2 -4c 2 (t-^) 2 + ^- 



w^-dTdt 
av 



'o Jr Q 
+ r(a - Ad - 2c) 



2r 2 |/f -4c^(t 



10 /j • udTdt 



where a = Ad — 2c — 1 + k for < k < 1 is a constant and E w is defined as follows: 



(2.11) 



E w (t) 



[w 2 (x,t) + w 2 (x,t) + \\7w(x,t)\ 2 }dtt 



An immediate corollary of the estimate is the following (Theorem 6.1 in |29j ) 

Corollary 2.2. Under the assumptions in Theorem \2. II) . the following one-parameter 
family of estimates hold true, for all r sufficiently large, and for any e > small: 



(2 



12) BT\ w + 2 [ [ e 2Tlfi F 2 dQ + const v [ [ F 2 dQ > k v [E w (0) + E W {T)} 
Jo Jn Jo isi 



for a constant k v > while BT\ W is given by: 



(2.13) BT\ W = BT\ W + const v 



dw 
dv 



w t 



dTdt 



h 

to J To 
+ 



w 2 dYr)dt 



0< 

T 

JTo 



ww t \dT dt 
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Remark 2.3. For the proof of the above Carleman estimate and the corollary, we refer 
to ESI and we omit the details here. 



2.2. Continuous Observability Inequality. Using the Carleman estimate in last 
section, we can prove the following observability inequality: 

Theorem 2.4. Under the main assumptions (A.l) and (A. 2), for the following initial 
boundary value problem 

' Wu(x, t) = Aw(x, t) + q(x)w(x, t) in Q x [0, T] 



(2.14) 



io(-, -|) = w (x) in f2 

w t (-, ?) = Wi(x) in f2 

frOM) = onr lX [0,T] 

IfKM) =g(x,t) onT x [0,T] 



where w G H^Sl), w 1 G L 2 (fi) ; # G L 2 (r x [0,T]) and q G L°°(0). VFe nave the 
following continuous observability inequality: 

II II 2 _i_ II II 2 <r H ( II II 2 _i_ II II 2 _i_ II II 2 I 

\\ w o\\m{n) + IPiIIi, 2 (q) ^ u ^ll w llL 2 (rox[o,T]) + ll w; tlli 2 (rox[o,T]) + llfl , llL 2 (r x[o,T]) J 

where T is as in (11.211) and C = C(fl, T, To, (p, t, q) is a positive constant. 

Proof. For the case when g = 0, we refer to [25] where the continuous observability 
inequality is established for zero Neumann data on the whole boundary. Here we give 
the proof for the case of general g G L 2 (T x [0, T]), which is still based on the proof 
in [29]. We first introduce the following result that is from the section 7.2 of |28j . 

Lemma 2.5. Let w be a solution of the equation 

(2.15) w tt (x, t) = Aw(x, t) + q(x)w(x, t) + f(x, t) in Q = Q x [0, T] 
with q G L°°(Q) and w in the following class: 

( w G L 2 (0, T; H 1 ^)) fl H 1 ^, T; L 2 (fi)) 

(2 ' 16) { ^,^GL 2 (0,T;L 2 (r)) 

Given e > 0, e > arbitrary, given T > 0, there exists a constant C — C(e, e , T) > 
such that 



T- 




( dw\ 



(2-17) / / |V,„ (1 «1^/IV/ < C ( / M + ^jj dr ^+l^lll 2( o,T ; ^-o ( n)) 

+ ll/ll^ +E o (Q) 

Now to prove (12.41) . we first establish the following weaker conclusion under the 
assumptions (A.l) and (A. 2) 

(2.18) e(-] < c( [ T [ [w 2 + w 2 + g 2 ]dT dt + \\w\\ 2 ^ ] 
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/ u 2 dVL < c\ 


1 \Vu\ 2 dQ + I u 2 dT 


; f u 2 dT <c 2 [ 


IV. 


Jn Jt 


Jfo Jn 



which is the desired inequality (12 .4p polluted by the interior lower order term \\w\\. 
To see this, we introduce a preliminary equivalence first. Let u G then the 

following inequality holds true: there exist positive constants < k\ < k 2 < oo, 
independent of u, such that 

(2.19) ki I [u 2 + \Vu\ 2 }dQ < [ \Vu\ 2 dQ + [ u 2 dT < k 2 [ [u 2 + |Vm| 2 ]^ 
Jn Jn Jf Jn 

where T is any (fixed) portion of the boundary T with positive measure. Inequality 
(I2.19P is obtained by combining the following two inequalities: 

(2.20) 

The inequality on the left of ( 12.201) replaces Poincare's inequality, while the inequality 
on the right of f )2.20p stems from (a conservative version of) trace theory. Thus, for 
w G H 2 (Q), if we introduce 

(2.21) e(t)= [ [\Vw(t)\ 2 + w 2 (t)]dQ+ [ w 2 {t)dT l 

Jn Jt 

where r = T \ Ti is as defined in the main assumptions, then (I2.19p yields the 
equivalence 

(2.22) aE(t) < e(t) < bE{t) 

for some positive constants a > 0, b > 0. 

Now in a standard way, we multiply equation (I2.15P by w t and integrate over Q. 
After an application of the first Green's identity, we have 

(2.23) ^ Qf [w t a + |V«;| 2 ]^ + jT w 2 dT ^j = jf ^w t dT + jT ww t dT 

+ / [q(x) + f]w t dtt 
Jn 

Notice that on both sides of f)2.23p we have added term f w 2 dT = f ww t dT . 

2 ^ J To Jr 
Recalling e(t) in (12.221) . we integrate (12. 23ft over (s,t) and obtain 



(2.24) e{t) = e(s)+2 



dw 
r dv 



w t dY + / ww t dTo 
r 



dr + 2 [ [ [q(x) + f]w t dQdr 
Js Jn 



We apply Cauthy-Schwartz inequality on [q(x) + f]w t , invoke the left hand side 

E(t) < -e(t) of (12391 . and obtain 
a 

(2.25) e(t) < [e(s) + N(T)\ + C T J e{r)dr 



(2.26) e(s) < [e{t) + N(T)] + C T / e(r)dr 
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where we have set 



2 ' 

e(|)+e(0) r t A7/ _ N 
-12^ K -^ e -°T a - jV(T) 



(2.27) JV(T)= / / fdQ + 2 [ [ ^w t dTdt + 2 [ [ \ww t \dT dt 

Jo Jn " Jo Jr ov J J Tq 

Gronwall's inequality applied on (I2.25p . (I2.26P then yields for < s < t < T, 

(2.28) e(t) < [e(s) + N(T)}e CT{t - s) ; e(s) < [e(t) + N(T)}e CT{t - s) 

We consider the following three cases here: 

Case 1: < s < i < ^. In this case we set t = -| and s = t in the first inequality 
of (12.281) ; and set s = in the second inequality of (12.281) . to obtain 

(2.29) e(|) < [e(t) + N{T)]e c ^- e(0) < [e(t) + iV(T)]e c ^ 
Summing up these two inequalities in (12.29!) yields for < t < 

<t) > 

(2-30) 2 

>^[^(2) + ^(°)] e ^ T¥ - Ar ( T ) 

after recalling the left hand side of the equivalence in (I2.22p . 

Case 2: ^ < s < t < T. In this case we set t = T and s = t in the first inequality 
of (I2.28P ; and set s = -| in the second inequality of (I2.28p . to obtain 

(2.31) e(T) < [e(t) + N(T)}e c ^; < [e(t) + N(T)}e c ^ 
Summing up these two inequalities in (12 .311) yields for ^ < t < T, 

e(t) > £ ® + < T K -o^ _ N{T ) 

(2.32) 2 

>pJ{Z) + E(T)]e-°*$-N(T) 

after recalling the left hand side of the equivalence in (I2.22P . 

Case 3: < s < ^ < t < T. In this case we set t = and s = t in the first 
inequality of (12.281) ; and set s = -| in the second inequality of (I2.28p . to obtain 

(2.33) £ (0) < [e(t) + N(T))e c ^; e(|) < [e(t) + N{T)]e c ^ 
Summing up these two inequalities in (I2.33P yields for ^ < t < T, 

s(t) > £ M±m e -c T z _ N{T) 
(2-34) 2 

>^[E(-) + E(0)]e- c ^-N(T) 

after recalling the left hand side of the equivalence in (I2.22p . 
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In summary, we get for any < t < T, 

(2.35) e(t) > \E{^)e- c ^ - N(T) 

We now apply the Corollary 2.2 of the Carleman estimate, except on the interval 
[e, T — e], rather than on [0, T] as in (I2.12p . Thus, we obtain since / = 0: 



(2.36) BT\ [e>T „ e]xT > k v E(e) 

where 5T|[ eT _ e ] xr is given as in (I2.13p . Since we have = on Ti x [0,T] and 
Qui 

-q— = g(x, t) on To x [0, T] by (I2.14p . with the additional information that h- v = on 
Ti by the assumption (A.l). Thus, by using the explicit expression ( 12.1 Op for BT\ W , 
we have that BT\u T - e ]xr is given by: 

(2.37) 

£T| [e , r -e]xr = 2r f [ e 2rip (w 2 - \Vw\ 2 )h ■ udTdt 

Je Jt 

+ 8cr I [ e 2rip (t - -)w t gdTdt 
Je Jf 2 

+ 4r f [ e 2rip (h-Vw)gdTdt 
Je Jr 

+ 4r 2 jf e 2 ^ (j/f - 4c 2 (t - T -) 2 + wgdTdt 



\t 2 (\h\ 2 - Ac 2 (t - ^) 2 ) + r{a - Ad - 2c) 



27 / I e 2TV 



w 2 h ■ vdTdt 



+ consty 



pi — e /• t>ti p r r 

/ / \gw t \dTdt+ / / w 2 dT dt + / / \ww t \dT dt 

_J e J Fq J J To J e J Fq 

Next, by the right side of equivalences ( I2.22p and ( I2.35p . we obtain 

(2.38) E(e) > ^ > (^j e~ c ^ - 2 jf J^\gw t \dTdt - 2 jf jf |wiu t |dr cft 

recalling JV(T) in ([2T2TD . We use (T2T38]) in f l2T36|) . Finally, we invoke estimate f[2~TTj) 
of Lemma 12.51 on the first and the third integral terms of ( 12. 37ft . This way, we 
readily obtain ( 12.181) . which is our desired inequality polluted by \\w\\ 2 x , 

To eliminate this interior lower order term, we can apply the standard compact- 
ness/uniqueness argument (e.g. [21]) by invoking the global uniqueness Theorem 7.1 
in ESI. □ 



inverse problem for a structural acoustic interaction 

3. Proof of Theorem 11.41 
We let w = w(f) = w t (f) then from (11 .2h we have w, u satisfy 
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(3.1) 



w tt (x, i) — Aw(x, t) — 


^(x)^ 


a;,*) 


= f(x)Rt(x,t) 


in ft x [0, T] 










on r x x [0,T] 


w(x,t) = —u tt (x,t) — 


A 2 u(x 


*)- 


A 2 u t (x, t) 


on T x [0, T] 


u(x,t) = $£(x,t) = 








on <9T x [0,T] 


< ^(x,t) = u tt (x,t) 








on T x [0,T] 


™(-,f) = 








in ft 


u>t(-,f) = f(x)R(x,^ 


) 






in ft 


«(-,?) = 








on r 


k «t(-,f) = 








on r 



Under the assumptions in Theorem 11.41 we can apply the Carleman estimate to the 
wave equation in the system (13.11) w tt (x,t) — Aw(x,t) — q(x)w(x,t) = f(x)R t (x,t) 
and get 

+ 2 [ e 2T v\fRt\ 2 dQ + C liT e 2T<7 [ w 2 dQ > (rep-2C T ) / e 2rip [w 2 t +\Vw\ 2 }dQ 
Jq Jq Jq 

+ [2t 3 /3 + 0{t 2 ) - 2C T ] [ e 2rip w 2 dxdt - c T T 3 e - 2T5 [E a (0) + E^T)] 

jQ(a) 

where the boundary terms are given explicitly by 

BT\u, = 2t f [ e 2rip (w 2 - \Vw\ 2 )h ■ vdTdt 

JO JTn 



'o Jr 

+ 8cr [ [ e 2rLp (t - ^)w t ^dTdt 
Jo Jr 



dv 



(3.2) 



+ 4r^ Je 2TV (h- Vw)^dTdt 
e 2r * (\h\ 2 -4c 2 (t- 



2 2r dv 



>o Jr 
2t f [ e 2 ^ 

Jo Jr 



2r 2 [\h\ 2 -Ac 2 (t-^) 2 



+ r(a -Ad- 2c) 



w 



2 h ■ udTdt 



Since we have the extra observation that Uu{x, t) = on To x [0, T] and note that the 
initial conditions u(x, ~) = Ut(x, -|) = on To, thus by the fundamental theorem of 
calculus we have u(x, t) = on T x [0, T] and hence from the coupling in the system 
(13~TT) we get 

(3.3) w(x, t) = -u tt (x, t) - A 2 u(x, t) - A 2 u t (x, t) = on T x [0, T] 
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and 

dw 

(3.4) — (x, t) = u tt {x, t) = on T x [0, T] 

Plugging (HOD and (E3D into (E2]), note also that — = on Y x X [0, T], therefore 
we get BT\ & = 0. 

In addition, in view of (ll.22p . (I1.23p . we have \fRt\ < C|/| for some positive 
constant C depend on R t . Moreover, notice that lim r 3 e~ 2rS = 0. Hence when r is 

sufficiently large, the above Carleman estimate can be rewritten as the following: 

(3.5) Ci, T / e 2rv5 [^ 2 + |Vw| 2 ]dQ + C2,r / e 2Tlp w 2 dxdt < C [ e 2rip \f\ 2 dQ + Ce 2T(7 

Jq ' Jq{ct) Jq 

where we set 

(3.6) Ci )T = rep - 2C T , C 2iT = 2r 3 /3 + £>(r 2 ) - 2C T 

and C denote generic constants which do not depend on r and henceforth we will use 
this notation for the rest of this paper. In addition, note that / is time-independent, 
so if we differentiate the system (13. ip in time twice, we can get the following wave 
equations for w t and w u : 

(3.7) (w t ) tt (x, t) - Aw t (x, t) - q(x)w t (x, t) = f(x)R tt (x, t) 
and 

(3.8) (w tt ) tt (x, t) - Aw tt (x, t) - q(x)w tt (x, t) = f(x)R m (x, t) 

Notice the assumptions (ll.22p . f 1 1 . 2 5 j) . therefore we have similarly as (13. 5 p the follow- 
ing estimates for the two new systems: 

(3.9) 

Ci, T / e 2rip \w 2 tt + \Vw t \ 2 ]dQ + C 2 , T [ e 2r *w 2 t dxdt <C I e 2rip \f\ 2 dQ + Ce 2ra 
Jq ' Jqm Jq 



and 
(3.10) 



Ci, T / e 2T ^[w 2 tt + \Vw tt \ 2 ]dQ + C 2)T [ e 2T ^w 2 t dxdt < C [ e 2Tlp \f\ 2 dQ + Ce 2ra 
Jq ' JqW) Jq 



where r is sufficiently large and Ci jT , C 2 , T are defined as in (13. 6p . 
Adding (13. 5p . (13. 9p and (I3.10p together we then have 

(3.11) C 1>T / e 2rip [w 2 t + w 2 tt + ™ttt + l v ^! 2 + l v ^*| 2 + \^w tt \ 2 ]dQ 
Jq 

+ C 2 , T [ e 2r *\w 2 + w\ + w 2 t ]dxdt <C[ [ e 2rLp \j\ 2 dQ + e 2ra ) 
Jq(u) \Jq J 



INVERSE PROBLEM FOR A STRUCTURAL ACOUSTIC INTERACTION 17 

Again we use the wave equation w tt (x,t) — Aw(x,t) — q(x)w(x,t) = f(x)R t (x,t), 
plugging in the initial time of t = ^ and use the zero initial conditions of w(-, ^) = 0, 
we have 

rj~\ rj~i rj~\ fjn rj~i 

(3.12) w u (x, — ) - AtD(x, — ) - q(x)w(x, — ) = w tt (a:, — ) = f(x)R t (x, — ) 

Since |i?t(x, ~)| > ri > from (11.231) . therefore we have |/(x)| < C|w it (x, ^)| and 

hence we have the following estimates on / e 2Tip \f\ 2 dQ: 

JQ 



(3.13) 

Jo Jn 



<C f T [ e 2T ^\w tt (x,^)\ 2 dndt 



o Jn 



Jn 1 



C {LJ 2 ^ e2Tl " (X ' S) l^( x ' s )| 2 )^ + J n e 2T ^ x ' 0) \w tt (x,0)\ 2 dn 
C Ucr jf jT a (| - S )e 2 ^ s V 4i (x, S )| 2 <M] 

T 

+ 2 / / 2 e 2 ^|^(x, S )||^(x, S )|d S ^+ / e 2T ^\w tt (x,0)\ 2 dn 
Jn Jo Jn 

T T 

<C\r I f" e 2 ^\w tt \ 2 dtdQ + f f" e 2 ^(\w tt \ 2 + \w m \) 2 dtdQ 
In Jo JnJo 

I \w tt (x,o)\ 2 dn 

n 

<C[t I e 2 ^\w tt \ 2 dQ+ I e 2 ^(\w tt \ 2 + \w m \ 2 )dQ 



= C((r+l) / e 2 ^\w tt \ 2 dQ + [ e 2 ^\w m \ 2 dQ 
V ./q Jq 

where in the above estimates we use the definition (12. ip and the property (12.41) of tp 
as well as Cauthy-Schwartz inequality. Collecting (I3.13P with (I3.1ip . we have 

(3.14) C 1>T / e 2T *[w 2 + w 2 t + w 2 tt + \Vw\ 2 + \Vw t \ 2 + \Vw tt \ 2 ]dQ 
Jq 

+C 2>T [ e 2T *{w 2 +w 2 +w 2 t ]dxdt <c((t+1) [ e 2T ^\w tt \ 2 dQ + [ e 2rLp \w m \ 2 dQ + e 
Jqm \ Jo Jq 
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Note that in (I3.14p . the right hand side term C J Q e 2Tip \wttt\ 2 dQ can be absorbed by 
the term Cx iT J Q e 2Tip [w 2 + w 2 t + w 2 tt ]dQ on the left hand side when r is large enough. 
In addition, since e 2rip < e 2ra on Q \ Q(a) by the definition of Q(cr), we have 

(3.15) 

C(r + 1) / e 2T *\w tt \ 2 dQ = C(r + 1) ( [ e 2T *\w tt \ 2 dtdx + / e 2 ^\w tt \ 2 dxdt 



< C(t + 1) ( f e 2rip \w tt \ 2 dtdx + e 2ra [ \w tt \ 2 dxdt 

\Jq(<t) JQ\Q(<t) 

<C(r + l) I e 2Tip \w tt \ 2 dtdx + C(T + l)e 

Jqm 



'Q(o-) JQ\Q{a) 
e ■ \wtt\ uiux" -|- o ^/ -tl)C 2 '" 

Again C(r + 1) Jq^ e 2TV \w t t\ 2 dtdx on the right hand side of (I3.15P can be absorbed 
by C% T $Qt a \ e 2Tip [w 2 + wf + w 2 t ]dxdt on the left hand side of (13.141) when taking r 
large enough. Therefore (13. 14j) becomes to 

(3.16) C' hT [ e 2 ^[w 2 + w 2 tt + w 2 m + \Vw\ 2 + |V^| 2 + \Vw tt \ 2 ]dQ 

JQ 

+ C' 2r [ e 2T *[w 2 + w] + w 2 tt ]dxdt < C ((r + l)e 2TCT + e 2rff + tV 2 " 5 ) 

jQ{a) 

Where we have 

(3.17) C' liT = rep-C, C'^ T = 2r 3 /3 + 0(r 2 ) 

Now we take r sufficiently large such that C 1 T > 0, C 2 T > 0. Then in (13.161) we can 
drop the first term on the left hand side and get 

C 2 T [ e 2rip [w 2 + w\ + w 2 t ]dxdt < C[(t + \)e 2ru + e 2ra \ 

(3.18) ' JQ{a) 

< C(t + 2)e 2ra 

Note again from (12 .6[) the definition of Q{<j), we have e 2Ttp > e 2l ~ a on Q(a), therefore 
(I3TT8D implies 

(3.19) C' 2 T I [w 2 + w\ + w 2 tt ]dxdt < C{t + 2) 

JQ{a) 

Divide r + 2 on both sides of ( 13.191) . we get 



(3.20) 



^- / [w 2 + w 2 + w 2 t ]dxdt < C 



By ( I3.17p . ^i. oo as r -> oo, thus (I3.20p implies that we must have w = on 
Q(cr) and hence we have 

(3.21) f(x)Rt(x,t) = w u (x,t) - Aw(x,t) - q(x)w(x,t) = 0, (x,t) e Q(a) 
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Recall again that \Rt(x, j)| > r\ > from (11. 23 j) and the property that Q D Q(cr) 3 
[t , ti] x n from fEZ77|> . Thus we have from (l3~2Tj) that f(x) = 0, for all x G ft. □ 

4. Proof of Theorem 11.51 

Setting f(x) = q(x) — p(x), w(x,t) = z(q)(x,t) — z(p)(x,t), u(x,t) = v(q)(x,t) — 
v(p)(x,t) and R(x,t) = z(p)(x,t), we then obtain (jl.2p after the subtraction of (II. ip 
with p from ( II. ip with q. Since R(x,^) = z(p)(a?, ?) = z (x) and R t (x,^) = 
Zt(p)(x,~) = zi(x), the conditions (I1.29P imply (ll.23p . In addition, the condition 
v(q)(x,t) = v(p)(x,t), x G r , t G [0,T] implies that u(x,t) = on r x [0, T] 
and (ll.30p implies ( 11.25jl . Therefore from the above Theorem 11.41 we conclude 
f(x) = q(x) — p(x) = 0, i.e., q(x) = p(x), iGSl. □ 

5. Proof of Theorem 11.61 

In relation with this system (13. ip . we define if) which satisfies the following equation 

>it(x, t) = Aip(x, t) + q(x)ip(x, t) in Q x [0, T\ 

on ri x [0,T] 
on T x [0,T] 
^(•,|)=0 in ft 

M;i) = fWR{xZ) in ft 

if), then we have y satisfies the following initial-boundary value problem 

f y u (x,t) - Ay(x,t) - q(x)y(x,t) = f(x)R t (x,t) in ft x [0,T] 

l(x,t) = onrx[0,T] 

= in ft 

= in ft 



(5.1) 

Set y = w 
(5.2) 







t) — Utt(x, t) 



9y_( 

dv 

y(-,|) 



It is easy to see that both (15 .ip and (I5.2p are well-posed. For the system (15.11) . we 
apply the continuous observability inequality in Theorem 12.41 to get 



(5-3) ||/i2(-,~)||i, (n) 



< C 



L2(r x[0,T]) 



+ 



^llL2(r x[0,T]) 



+ II ^ I 
dv 



li 2 (r x[o,T]) 



Notice that \R{x, f )| > r > 0, §£(x, t) = u tt (x,t) on T x [0,T] and §J(:M) = on 
r x x [0,T], therefore we have from (15. 3p 

(5-4) H/Hi^n) < C (||^||L2(r x[o,T]) + ||^IU 2 (r x[o,T]) + IKt IU 2 (r„x[o,T])) 
On the other hand, for the system (I5.2p . we have the following lemma: 

Lemma 5.1. Let q G L°°(ft) and R(x,t) satisfies R t G H^+ e (0,T; L°°(ft)) for some 
< e < \ as in Theorem \l.b\ If we define the operators K and K\ by K, K\ : 
L 2 (ft) -> L 2 (r x [0,T]) ; suc/i toot 

(5.5) (Kf){x,t) = y(x,t), {K x f){x,t) = y t {x,t), xeT ,te [0,T] 

where y is the unique solution of the equation (15. 2p . Tnen if and i^i are froin compact 
operators. 



(5.9) 
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Proof. It suffices to just show that K\ is compact, then it follows similarly that K is 
also compact. Since / G L 2 {VL) and R t G #2+ e (0,T; L°°(ft)), we have 

(5.6) fR t EH^(0,T;L 2 (n)) 
Therefore we have the solution y satisfies (e.g. Corollary 5.3 in [27J) 

(5.7) !/eC([0,n^)), y t eC([0,T};H^(Q)) 
Hence by (15. 6p . q G L°°(Q) and y tt = Ay + q(x)y + fR t we can get 

(5.8) y a eL\0,T;H-^(Q)) 

In addition, by ( 15 .7|) and trace theorem we have y t G C([0, T]; H e (F)). Since the em- 
bedding H e (T) — > L 2 (T) is compact, we have by Lions-Aubin's compactness criterion 
(e.g. Proposition III. 1.3 in [33]) that the operator K\ is a compact operator. □ 

Now we have that the inequality (15.41) becomes to 

L 2 (Q) < C (||V , IU 2 (r x[o,T]) + ||^t|U 2 (r x[0,T]) + \\Utt\\L 2 (r„x[0,T])) 

< C (\\w - t/||L2(r x[o,T]) + \\w t - yt\\L 2 {r x[o,T]) + IKt|U 2 (r ( ,x[o,T])) 

< C {\\w\\l 2 (t x[q,t]) + ||^||L 2 (r x[o,T]) + H M «IU 2 (r x[o,r])) 
+ C||?/||L 2 (rox[o,T]) + C||2/jL 2 (r x[o,T]) 

= C (||^||L2(r x[o,T]) + ||^IU 2 (r x[o,T]) + ||wtt|U 2 (r x[o,T])) 

+ C\\Kf\\ L 2 (r x[0,T]) L 2 (r x[0,T]) 

< C (||wtt||L 2 (r x[o,T]) + ||Mttt||L 2 (r x[o,r]) + ||^ 2 w«||L 2 (r x[o,r])) 

+ C||^/||L 2 (r x[0,T]) + C||^l/||L 2 (r x[0,T]) 

where in the last step we use the coupling w(x, t) = —u u (x, t) — A 2 u(x, t) — A 2 u t (x, t) 
on T x [0,T] from (13.11) and again the initial conditions «(•, f-) = Ut(-, ?) = on 
T x [0, T] so that by the fundamental theorem of calculus, we have 

(5-10) |MU 2 (r x[o,T]) < C|ht||L 2 (r x[o,T]) < C\\u tt ||L 2 (r x[o,T]) 

To complete the proof, we need to absorb the last two terms in (15. 9p . To achieve 
that, we apply the compactness- uniqueness argument. For simplicity we denote 

\\ u \\x = ||^tilU 2 (r x[o,T]) + IKtt||L 2 (r x[o,T]) + II A 2 Wf t || L 2( roX [ 0iT ]) 

Suppose contrarily that the inequality ( 1 1 . 3 3 [) does not hold. Then there exists /„ G 
L 2 (f2), n > 1 such that 

(5.11) ll/n|U 2 (fi) = h n>l 
and 

(5.12) lim |K/„)|U = 

From (15. lip , there exists a subsequence, denoted again by {/„} n >i such that /„ 
converges to some fo G L 2 (Q) weakly in L 2 (Q). Moreover, since K and K\ are 
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compact, we have 

(5.13) lim \\Kf n - Kf m \\ L 2 {r x[QT]) = 0, lim \\K x f n - #i/ m ||i,a(r x[o,r|) = 

m,n— >-oo m,ra— >oo 

On the other hand, it follows from (15.91) that 

Wfn fm |U 2 (f}) < C\\u(f n ) - u(f m )\\ X + C\\Kf n - i^/m||L 2 (r x[0,T]) 

, , + C\\K l f n - Kif m \\ L 2(r x[o,T\) 

( ' K 1 <c|| w (/„)|U + c'||«(/j|U + c'||x/ n -x/ m || i2(roX[ o, T]) 

+ C\\Kif n — i^i/ TO ||L2(r x[Q,T]) 
Thus by (I5.12p and (15. 13)) . we have that 

(5.15) lim ||/ n - / m ||£2( ) = 

m,n— >oo 

and hence f n converges strongly to fo in L 2 (Q). So by (15.111) we obtain 

(5.16) ||/o|b(n) = 1 

On the other hand, by (11.231) and a usual a-priori estimate, we have that 
^ \\w(f)\\c([o,T];m(n)) + \\wt(f)\\c{[o,T];Li{n)) < C\\fRt\\ L i(o }T . L 2^)) 

< C||^t||Ll(0,T;X«>(n))||/||L 2 (n) 

Hence trace theorem implies that 

(5-18) ||tS(/)|| tfC r x|p l r|)<C7||/|| L3(O ) 

where C > depends on \\R t ||Li(o,T;L°°(n))- Therefore by (15.181) we have 

(5.19) lim ||«3(/ ft ) - w(/o)||L2(r x[o,r]) < C lim ||/ n - / ||x2(n) = 

n— >oo n- »oo 

Moreover, by (15.121) and the coupling w(x, t) = —u a (x, t) — A 2 u(x, t) — A 2 u t (x, t) on 
r x [0, T], we have 

(5.20) lim ||w(/„)|| L 2 (rox[OT]) < lim \\u\\ x = 

n— >oo n— yoo 

Thus by (15TT9D and (Ex20D . we obtain 

(5.21) w(f )(x,t) = 0, xeT ,te[0,T] 

Therefore from (13.11) we have u = u(fo) satisfies the initial boundary problem: 

' -u tt (x, t) - A 2 u(x, t) - A 2 u t (x, t) = in T x [0, T] 

u(x,t) = % (x,t) = onar x[0,T] 

w(-,~) = in T 

k M t (-,|) = in T 

which has only zero solution, namely, we have u(fo)(x,t) = 0,x G T ,t e [0,T]. 
Therefore by the uniqueness theorem 11.41 we have fo = in which contradicts 
with (I5.16p . Thus we must have 

(5-23) < C [\\uu\\l2(Tox[o,t]) + ||^i||L2(r x[o,T]) + || A 2 u tt \\ L 2 {rQX[0tT]) ) 

and the proof of the theorem is complete. □ 



(5.22) 
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6. Proof of Theorem 11.71 
We now go back to the original system (11. ip . 

Case 1: n=2. Let z = z t and v = v t , then the system (11.11) becomes to 



(6.1) 



dv 

dv 



in tt x [0, T] 
on ri x [0,7] 

v a (x, t) - A 2 v(x, t) - A 2 v t (x, t) on T x [0, T] 
= 

(x 
zi{x) 



z a {x, t) = Az(x, t) + q(x)z(x, t) 
Zt(a;,t) 

dz 



|) = A;z (x) + g(x)z 



t u t (-, |) = -zi(ar) - A 2 w (x) - AVl(x) 



on ar x [0,T] 
on T x [0,T] 
in f2 
in Q 
on r 
on r 



By using the similar operator setting as in section 1.2 and notice the new initial 
conditions, we can compute the domain of the operator A 2 '. 

(6.2) 

D(A 2 ) = {[zq, zi,« , v x ) T : (zi, (-A N + q)z + Bv u v u -A(v + Vx) - £*Zi G £>(-4)} 
= {[z ,Zl,i>o,Wif : Zi G fl? 1 (n),(-^ JV + g),2b + -BT;i G ^(fi),^ G # 2 (r c 
«! G iJ 2 (r ),A^ + v x ) + 5*^ G H 2 (T ), 



Oh 



v 1 -A(v + v 1 )-B*z 1 eD(A),^ 



r = -A(«b + «i) - 
•9^0, 



{[^0,^1, v ,^i] T : z\ e H^(Q), (A + g)z e H^(Q), -^j|r = t>i, 

« e // 2 (r )^i g # 2 (r ), A(^ + Wl ) + g # 2 (r ), 



«i - A{v + vx) - B*zi G D(A) 



9^1 



I To 



-A(v + Vi) - B*2 X } 



{[zb,zi,i>b,vif : e H^{Q),Zx G H^(tt),Vo G # 2 (r ),vi g # 2 (r ), 
A( Vo + m) + B* Zl g # 2 (r ), Vl - A(^ + vx) - B*z x G £>(A), 



9^ 



I To 



Vx, 



dzx 
dv 



r = -A(v + Vx)-B*zx} 



where in the last step zo G Fp (fi) is from elliptic theory when provided that g(x) G 
WA°°(ft). Therefore when z G H^(Q), z x G #f\(ft), t>o G #o( r o), «i G # 2 (r ) with 
compatible conditions as in L>(.4 2 ) and g G PV 1 ' 00 ^), then from semigroup theory 
we have that the solution of (16. ip satisfies 

(6.3) z t G C([0,T\;H\n)), z tt G C([0,7]; L 2 (Q)) 
Hence we have on the one hand 

(6.4) ^il^T;^)) 
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On the other hand, from (I6.3P and z u (x, t) = Az(x, t) + q(x)z(x, t), we have by elliptic 
theory that 

(6.5) zt = z G L 2 (0,T;H 2 (n)) 
Interpolate between (16 .41) and (16. 5ft . we have for < e < |, 

(6.6) zt E H^ +e (0, T; Hi- e (Q)) C F^ +£ (0, T; L°°(fi)) 

where the inclusion is by Sobolev embedding theorem. 

Case 2: n=3. We let z = z tt , v = Vu, then we have z, v satisfy 



zu(x, t) = Az(x, t) + q(x)z(x, t) 



(6.7) 



dz 



M t (x,t) - 
v(x, t) - 







-v tt (x, t) - A 2 v{x, t) - A 2 v t {x, t) 
x,t) = 



du 



%{x,t) = v t {x,t) 



du 

T 
2 

2 > 
T\ 



US 



Az (x) + q(x)z 
-- Azi(x) + q(x)z 1 

—z\{x) — A 2 v {x) — A 2 vi(x) 
- —Az (x) — q(x)z (x) — A 2 vi(x) 

+A 2 z 1 (x) + A\(x) + A%(x) 



in Q x [0,T] 
on r x x [0,T] 
on T x [0,T] 
on dV x [0,T] 
on T x [0,T] 
in Q 
in Q 
on r 

on r 



Then still using the similarly operator setting as before we can compute the domain 
of^ 3 : 

(6.8) 

D(A 3 ) = {[z ,Zx,v ,vx] T : {z x ,{ r A N + q)zQ + Bv u v 1 ,-k{v Q + v 1 ) -B*z x e D(A 2 )} 
= {[zo^^vo^vxf : Zl G H^(n),(A + q)z G H^(Q),^\r = v h 
v G H 2 (r ), Vl G H 2 (T ),A(v + v 1 ) + B*z 1 G H 2 (T ), 
A( Vl - A(u + «i) - S*^) + B*[(-A N + q)z + B Vl ] G 5 2 (r ), 
A(w + «i) + B* Zl + A[ui - A(v + vi) - + B*[(-A N + q)z + Bv x ) G D{k) 



^|r = -A(u Q + ui) - 

OP 



d[{-A N + q)z + Bv x ] 
du 



|r 



- A[«i - A(v + Vl ) - B* Zl ] - B*[{-A N + q)z + Bvx}} 

= {[zo^x^o^xf : z G H^(n), Zl G ^(fi),^ G if£(r ),T>i G # 2 (r ), 

A(^ + Wl ) + B* Zl G H 2 ( r o), -^j|r„ = vx, — |r = ~A{v + v t ) - B*z x on r , 
A(v + «i) + + A[«i - A(u + vx) - B* Zl ] + S*[(-Ajv + q)z + G D(A) 
A(vx - A(v + vx) - B* Zl ) + B*[(-A N + q)z + Bv x ] G H 2 (T ), 

dK-A^q^ + Bvx] ^ = _ A(vo + v0 _ _ n{ _ AN + ^ + ^ ]} 
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where in the last step z G Hj, (Q) is from trace theory of solving ^|r — v± E 

H 2 (T ). Therefore when z G H^(n), z x G H^(Q), v G H$(T ), ^i e #o( r o) with 
compatible conditions as in D(A 3 ) and q G H/ 2 '°°(fi), then from semigroup theory 
we have that the solution of (16. ip satisfies 

(6.9) % G C([0, T]; H\Q)), =z tt G C([0, T]; L 2 (fi)) 
Hence we have on the one hand 

(6.10) z t G H 2 (0,T; H 1 ^)) 

On the other hand, from (16. 9p and z a (x, t) = Az(x, t) + q(x)z(x, t), we have by elliptic 
theory that 

(6.11) z tt = zeL 2 (0,T;H 2 (n)) 
which implies 

(6.12) z t G i^ 1 (0,T; H 2 (Vl)) 

Now interpolate between (I6.10p and (16.121) . we have for < e < |, 

(6.13) z t G Ft (0, T; /ft (O)) C #5 +e (0, T; L°°(0)) 

where the inclusion is again by Sobolev embedding. 

Hence in either case n = 2 or n = 3, under the assumptions on the initial data 
[zq, Zi, v , vi] and q(x), p(x) in Theorem 11.71 we have that z t G if2+ e (0, T; L°°(Q)). 
Thus when we again set /(x) = q(x) —p(x), w(x, t) = z(q)(x, t) — z(p)(x, t), u(x, t) = 
v(q)(x,t) — v(p)(x,t) and R(x,t) = z(p)(x,t) as in section 4, we obtain (ll.32p that 
Rt G -£f2+ e (0,T; L°°(Q)) and hence all the assumptions in theorem 11.61 are satisfied. 
Therefore, we get the desired stability (ll.34p from the stability f l 1 . 3 3 [) of the linear 
inverse problem in Theorem 11.61 □ 

7. Concluding remark 

As we mentioned at the beginning and the calculations of D(A 2 ) and D(A 3 ) show, 
the lack of compactness of the resolvent limits the space regularity of the solutions 
for the wave equation parts since we always have the elliptic problem for z or z t such 
that (A + q)z G L 2 {Vt) with g G #o( r o) provided q in some suitable space. Therefore 
the best space regularity that z could get is 2 + | = | from elliptic and trace theory. 
As a result, our argument of the stability in the nonlinear inverse problem will only 
work for dimension up to n = 7 as we need the Sobolev embedding Hz(Q) C L°°(Q) 
in order to achieve the space regularity of z t in L°°(Q) which is needed in the proof. 
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